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The field of topological insulators (TIs) is rapidly growing. Concerning possible applications, the 
search for materials with an easily controllable TI phase is a key issue. The quantum spin Hall 
effect, characterized by a single pair of helical edge modes protected by time-reversal symmetry, has 
been demonstrated in HgTe-based quantum wells (QWs) with an inverted bandgap. We analyze 
the topological properties of a generically coupled HgTe-based double QW (DQW) and show how 
in such a system a TI phase can be driven by an inter-layer bias voltage, even when the individual 
layers are non-inverted. We argue, that this system allows for similar (layer-)pseudospin based 
physics as in bilayer graphene but with the crucial absence of a valley degeneracy. 

PACS numbers: 73.20.-r, 73.43.-f, 73.21.Fg, 73.61.-r 
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Since the understanding of the topological nature of 
the quantum Hall efFectii^, topological phases have be- 
come one of the most active research fields in condensed 
matter physics. More recently, a new topological phase 
preserving time-reversal symmetry (TRS), the quantum 
spin Hall (QSH) phase'^'* has been discovered. The QSH 
phase has been theoretically predicted^' and experimen- 
tally realized in 2D HgTe/CdTe QWs^. The crucial in- 
gredient of this narrow gap semiconductor material is 
strong spin-orbit coupling, which determines the inverted 
band structure of HgTe. The experimentally accessible 
parameter tuning the band structure from normal (CdTe- 
like) to inverted is the thickness of the HgTe QW. One 
year later, 3D TIs supporting chiral fermions as surface 
states have been proposed and observecl2r~. These two 
phenomena are examples of the general concept of a TI 
which is a TRS preserving system with a bulk insulating 
gap which features topologically protected edge states 
due to the Atiyah-Singer index theorem which is in this 
context referred to as the bulk boundary correspondence. 
TI phases are characterized by a Z2 topological invari- 
ant'^ 



f3,13 



In this work, we extend the Bernevig-Hughes-Zhang 
(BHZ) model^ to account for a double QW (DQW) of 
HgTe (see Fig. 1). We analyze how the topological fea- 
tures of the DQW depend on a generic tunneling Hamil- 
tonian effectively connecting the two wells and on the 
applied voltage V between them. In particular, we show 
that a QSH phase can be driven by V also when the two 
QWs are individually trivial. We also derive a reduced 
2-band model that well captures the topological features 
of the system. We give analytical expressions of the ef- 
fective parameters that map the reduced model to the 
BHZ-model of a single HgTe QW [in Eq. (10)] as a func- 
tion of the tunneling matrix elements and voltage V . We 
show that a non-trivial Z2 topological invariant is accom- 
panied with the appearance of a single pair of helical edge 
states for which we calculate the energy dispersion and 
spinors. In Fig. [^fa) we show the phase diagram of the 
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FIG. 1: (Color online) (a) Schematic representation of the 
HgTe/CdTe DQW made by a front and back well with a finite 
overlap of their individual envelope functions Xo(-^) ^^'^ Xo(-2) 
as shown in (b). (c) Process of band inversion in a DQW 
caused by a potential bias. 



model as a function of tunneling amplitude and applied 
voltage. 



I. DOUBLE QUANTUM WELL MODEL 

The system under investigation consists of a DQW of 
HgTe as shown in Fig. [TJ which can be thought of as a 2D 
TI bilayer. The system is contacted to gates that control 
the potential bias V between the front (f) and back (b) 
QWs. The spectrum of a single HgTe-based QW near 
the F-point is effectively described by the BHZ model^ 



h{k) = 
d ^ 



( h{k) 
V h*{~k) 

d- a (1) 

(C - Dk^, Ak^, -Aky, M - Bk^) , 



where a are the PauH matricesi^ associated with the 
band-pseudospin degree of freedom (band Ex or ii\). 
Ho is represented in the basis {|i<^i+), |-ffi+), |-Ei— ), 
\Hi~)}, where the Ei states {J^ — ±1/2) are a mixing 
of the s-hke Fg band with the Fg hght-hole band, while 
Hi {Jz = ±3/2) is basically the Fg heavy-hole bandit. In 
the following we will use the parameter values A = 375 
meV nm, 5 = -1.120 eV nm^ and Z? = -730 meV nm^, 
estimated by a comparison with the 8x8 Kane Hamil- 
toniaii^, and assume C = without loss of generality. 
The Dirac rest mass M depends on the QW thickness 
and M < corresponds to the inverted (QSH) regime 
whereas M > corresponds to the normal regime. In 
first approximation, Hq is block diagonal in the Kramer's 
partner or spin degree of freedom^ (with s the vector of 
Pauli matrices). Because we consider only systems with 
TRS, we restrict ourselves to the block h{k), from which 
the results can be extended to the other one by applying 
the time reversal operator T — iSyK. 

The Hamiltonian h{k) in Eq. [T] describes the in-plane 
electronic motion inside a QW layer. Along the confine- 
ment direction of the QW (Z), electrons are described by 
the envelope function Xo{z) (integrated out in the BHZ 
model). When the two QWs of a DQW are sufficiently 
separated, an electron will be localized either on the f or 
on the b layer, with envelope function Xoi^) ^-nd Xo{z)i 
resp. When the two layers are placed close to each other, 
Xo{z) s-nd Xo(-^) acquire a finite overlap, accounted for 
by a tunneling Hamiltonian Ht- This description should 
take into account that xo(-z) is a spinor with components 
on the El and Hi bands. To first order in k, the tunnel- 
ing Hamiltonian has the form 



Ht 



+2(^p. 



3(A) • a 



2 2 

A = {Ao,akx,-aky,Az) , 



Vy 



(2) 



where Vi are Pauli matrices associated with the layer 
projection [layer pseudospin (LPS) V]. Ht generates 
bonding/antibonding states with energy splitting Aei — 
Ao + Az {Am = Ao - A^)^ of the El (HI) band. The 
lowest order off-diagonal term is ak± , due to the axial 
symmetry and the Jz character of El and HI bands- . 
In the Appendix A, we present a realistic estimate of 
HgTe DQW tunneling parameters. 

The DQW Hamiltonian for a single Kramer's block is 
therefore 



H = d.,Vo + ^^V.-^^:^Vy-iSd-a)Vz, (3) 



where Sd is due to the possible variation of the parame- 
ters of h{k) between the f and the b layer. In particular, 
we consider Sd = {V, 0, 0, 0), with V the interlayer bias. 



II. BAND STRUCTURE AND TOPOLOGY 

The topological properties of a fully gapped 2D TRS 
preserving system are described by the Z2-invariant v —. 
When the model Hamiltonian is block diagonal with re- 
spect to the Kramers partner spin s, the system can be 
thought of, from a topological point of view, as two copies 
of an anomalous quantum Hall effect, related by time 
reversali^. Each of the blocks is topologically character- 
ized by its first Chern number C^ and C^ resp., which we 
call the Kramers Chern numbers (KCNs). TRS immedi- 
ately implies that the two KCNs obey a zero sum rulei^, 
meaning C-f = — C^. The Z2 invariant is then given byi^ 

^ = ^^ ^ ^^ (mod 2) = C^(mod 2), (4) 

with the KCN defined as 



(5) 



2tt Jt-2 
and where the Berry curvature is given by 

Hk)= 5](rf«|)A(dK)), 



(6) 



where A stands for the exterior product. The sum is 
over the occupied bands and |u^) are Bloch states of the 
Kramers block Hamiltonian. In order to make the Chern 
number well defined in local models, like Eq. [31 a lattice 
rcgularization has to be included to compactify the k- 
space. In our case, since the curvature decays rapidly 
away from the F-point, the integral in Eq. ([5]) over the k- 
space M? converges stably towards the KCN of the lattice 
regularized model. 

Let us first analyze the bulk dispersion curves de- 
scribed by Eq. [3] to get an intuition for possible local 
topological phase transitions^ at the F-point, induced 
by band crossings. For fc = we obtain the four eigenen- 
ergies 



1 



^±,^(0) = TiM ± -V (Ao + TiAzY + 1/2, 



(7) 



with ?7 = ±1. Due to the interlayer tunneling, the El 
bands (centered ai E = M) suffer an energy splitting, 
at the F-point, of •y/AlJj+V^, while for the HI bands 
(centered at ^ = -M) the splitting is \/ A]j^ -I- V^. We 
define a local energy gap (LEG) at fc = as 



Eg{0) = 2M- 



sgn(M) 



i^^li+V' + ^J^li+V')- (8) 



When the LEG changes sign an inversion of bands with 
El and HI character occurs at the F-point, accompanied 
by a possible topological phase transition. For a fixed 
set of parameters, the condition Eg{0) = defines (if 
existent) the critical value of the bias Vc driving the phase 
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system in Fig. [2][b). These states have the pecuharity 
of supporting at the same time both oscihatory and de- 
caying behavior [clearly evident in Fig[3{a)], due to the 
contribution of modes with complex ky, which is related 
to the Mexican hat shape of the bulk bands. Approach- 
ing the bulk bands, the oscillations tend to become less 
pronounced, as shown in Fig[3Jb). The edge modes have 
the major contribution from the Hi band of the f layer, 
which consists in 60% of the probability weight near the 
conduction band and rises to 90% near the valence band. 
Under the reversal of V, the probability weights of the b 
and f layer interchange. 



FIG. 2: (Color online) (a) KCN of a DQW as a function 
of V and A^ with a 7^ and the physically plausible condi- 
tion Ao — Az, corresponding to Ahi ~ (see Appendix IA 2|l . 
The dashed line, indicating the threshold bias is obtained with 
Eq. |9l The internal region is topologically trivial with Cf = 
{Cf = 2) for M > {M < 0), while the outer one is topolog- 
ically non-trivial with Ci^ — 1. (b) and (c): first conduction 
and valence bulk bands and edge states (if present) disper- 
sion curves obtained by numerically solving Eq. [3] (full lines) 
and analytically from the reduced model (dashed line) with 
M = 6 meV and taking into account a spacing distance be- 
tween the two QWs of 6 nm, leading to the estimated value 
of A_Ei = 6 meV, a — 5 meV nm (Table |l] in the Appendix). 
The two situations correspond to the cross and the circle in 
(a), resp. For 1/ = 15 meV, the DQW is topologically non- 
trivial (b), while for V — 9 meV, the system is trivial with no 
edge states (c). 



III. REDUCED HAMILTONIAN 

In order to get an analytical understanding of the sys- 
tem, we derive a reduced 2-band low-energy Hamiltonian 
hii, with the assumption that the relevant energies are 
small compared to -j « ±M. hu is represented in the 
basis {\Ei,b),\HiJ)} or {\EiJ), \Hi,b)} for the + and 
— sign, resp. Similar to BLG, we separate the DQW 
Hamiltonian in Eq. [3] in a low-energy hu and a high- 
energy hfih part connected by the off-diagonal blocks 
hih = hhi- The reduced 2-band model is obtained then 
as hu = hu + huiGlf^hhi, where Gl,^ ^ {E ~ hhh)~'^ — ■ 
The reduced model hu can be mapped to a TI Hamilto- 
nian of a single layer (Eq. [1} , but with the renormalized 
parameters 



ransition: 
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(9) 

This is confirmed by the explicit calculation of the Z2- 
invariant v with Eq. |4l as shown in Fig. [2l where we 
calculated the KCN C^, as a function of A^ and V. C^ 
is insensitive to the value of a, as long as a and/or Az 
remains finite, in which case an insulating gap arises. 
Interestingly, for \V\ = |T4|, lowest bulk bands have at 
small k a Dirac-like dispersion, while a tunable gap is 
developed for \V\ < \Vc\. In the region \V\ < |Vc|, the 
system supports an even number of pairs of helical edge 
modes N = \C^\, in particular iV = for M > and 
A^ = 2 for Af < 0. In the outer region, instead, the 
DQW is topologically non-trivial and allows for a pair of 
time-reversal protected helical modes. An example of the 
DQW spectrum in the two regimes, described by Eq.[3l is 
shown in Fig[21Jb) and (c) (full lines) for a positive mass 
M = 6 meV and a bias of 1^ = 15 meV and 9 meV, resp. 
Edge states are obtained by numerically solving Eq. [3] 
with open boundary conditions. 

In Fig[3l^a) and (b), we plot the wave function of the he- 
lical edge states for E = 0.21 meV (near the Dirac point) 
and E = —1 meV (near the valence bands), resp. for the 
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while D = D remains unchanged. The bulk disper- 
sion is characterized by a Mexican hat dispersion for 



(M < ^^, where ^ 



sgn 



B{D^ - B^ 



with a gap 



E„ 



-^x/(^^ 



AMB){D^ - i?2). If the condition 
for the existence of the Mexican hat is not fulfilled, the 
dispersion has a minimum at /c = and a gap of 2M. For 
a barrier thickness of 6 nm, we obtain Eg ~ —1.6 meV 
(see Table im in the Appendix for details), which does not 
strongly depend on M. This suggests that the parame- 
ters obtained are in the correct ballpark for experiments. 
The bulk dispersion curves are quite well reproduced 
by hu, as demonstrated in Figs. [2jb) and (c) (black 
dashed lines). However, hu predicts edge modes with 

purely linear dispersions E = —~M + sAk^J ^~^, 

where s = ±1 denotes the two spin-blocks [colored 
dashed lines in Figs.[51Jb)]. Instead, edge dispersions nu- 
merically obtained from the full model [colored full lines 
in Figs, mb)] show a marked non-linear behavior when 
approaching the bulk bands, with an accompanying shift 
of the Dirac point. 
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FIG. 3: (Color online) Edge mode wave function along the 
Y direction, calculated from the full model with vanish- 
ing boundary condition at y — 0, for the DQW system of 
Fig. ^h) near the Dirac point E ~ 0.21 meV in (a) and for 
E ~ —1 meV (near the bottom of the valence bands) in (b). 
In particular, we show the behavior as a function of y of the 
projection on the Ei and Hi bands of the front (f ) and back 
(b) QWs. 



For the reduced model the KCN can be calculated an- 
alytically since J-" assumes the simple form^ 



-^/.ii/ = -yd ( df^d X d^d 



(11) 



where d is the unit vector of {Ak^, ~Aky, M — Bk^). The 
topological features of hu are therefore analogous to those 
of a single layer of 2D TI. However, the bias voltage en- 
ters the Dirac mass parameter M thus providing an ad- 
ditional knob to tune a system of two trivial insulators 
into a single- valley Dirac system {M = 0) and eventually 
into a nontrivial TI phase. Since the phase transition at 
M = is local at the F-point^ the analytical topological 
predictions of the reduced model are inherited by the full 
model and survive a lattice regularization which makes 
the KCN mathematically well defined. 

Finally, we note that the proposed mechanism of a 
voltage bias induced band inversion in a DQW is by no 
means restricted to HgTe based QWs. The method in 
principle even applies to large band gap semiconductors 
like GaAs (F ~ 1 eV) but the presumably large inversion 
asymmetry, not captured in our model, would have to be 
taken into account explicitly. 



IV. DISCUSSION AND CONCLUSION 

It has been proposed earlier that, in single QWs, the 
topological phase transition may be induced by an effec- 
tive applied potential^"-2i. In Refs. |22|24| . calculations 



show that in a single Hgi_2;Cd3;Te QW, a potential bias 
along the well thickness of the order of hundreds of meV 
[much larger than V needed in our proposal (on the order 
of a few meV)] induces a band inversion. The experimen- 
tal feasibility of this approach is arduous due to the large 
field required. In Ref. l23|, a type-II quantum well made 
of InAs/GaSb/AlSb, which has an intrinsically inverted 
bandstructure, is considered. The authors show that the 
spatial separation of El and HI bands offers a key to 
drive the transition from the TI to the normal phase with 
an applied gate bias. Recent experiments^i2S. provided 
the first evidences pointing towards the presence of a TI 
phase in these kind of structures. The proposal in Ref^ 
requires the growth of a type-II QW, which is a strongly 
asymmetric structure interfacing two specifically chosen 
semiconducting materials with different electron affinities 
(InAs and GaSb). On the contrary, we envision a DQW 
with each well made from the same material that might 
be more feasible to experimentally realize clean samples 
and offers the possibility of tuning the tunneling param- 
eters by modifying the barrier thickness or its chemical 
composition. This tunability could also be used to cre- 
ate contacts between two helical Luttinger liquids when 
each QW is in the non-trivial regime and V ~ 0. This 
contact can be either of electrostatic nature which allows 
to study Coulomb drag?^ or a new type of inter-edge 
correlated liquid^-, or can induce the tunneling between 
helical edge states^E. 

Further, we note that the LPS (V) [sec Eq. (2)], that 
operates on the wave function amplitude on the (/) and 
(6) layers, can be manipulated by the various parame- 
ters, in particular by the voltage V. As already men- 
tioned, in the regime V/2 « ±M, the low-energy physics 
of the DQW is described by Eq. (1), with parameters 
from Eq. (10), where the basis is given by eigenstates of 
the LPS with eigenvalues / and b resp. If M in Eq. (10) 
is zero, V points in-plane for small in-plane /c-vector and 
has a Berry phase of n like in graphene. A finite mass 
M tilts V out-of plane and starts to localize V in one 
of the layers. This physics is actually reminiscent of bi- 
layer graphene (BLG) in the presence of an interlayer 
voltage—, except that the low-energy physics is governed 
by a LPS with Berry phase 27r. We expect that inter- 
esting proposals that utilize the LPS in BLG, e.g. in 
a LPS valvals, also apply here, but with the important 
absence of a valley degeneracy that is potentially harm- 
ful for operating such devices^. Also, the application 
of a voltage domain in BLG predicts valley-filtered edge 
states at the domain- wallSi which in our case would result 
in (Kramers) spin-filtered edge states without valley de- 
generacy at a mass(M) domain — i.e. helical edge states. 
Such mass domains would allow the controlled creation 
of helical edge states not only at the physical sample 
boundaries of DQW systems but also in their bulk. 

In summary, we have investigated DQW structures us- 
ing the BHZ model of HgTe QWs and a generic tunneling 
Hamiltonian connecting the layers. An interlayer poten- 
tial bias on the order of the layer bandgap can drive a 



topological phase transition even if the individual QWs 
are in the normal regime. We calculate the Z2 topological 
invariant and the helical edge states which in a reduced 
model obtain simple analytical structures. These results 
suggest DQWs as potential candidates for an all tunable 
topological insulator or Dirac system which would have 
desirable properties for applications. 



where fn;k« {z) are the envelope functions, ky — [k^, ky) is 
the wave vector in the plane of the QW, and Unir) is the 
usual basis set for the eight-band Kane model^""— which 
is assumed to be the same in HgTe- and Hgi-^rCd^^Tc- 
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Appendix A: Estimate of the tunneling parameters 
in a HgTe DQW 

In the present Appendix, we will provide a quantita- 
tive calculation of the tunneling matrix elements [Eql5] 
between the electronic states belonging to the front 
(f) and to the back (b) HgTe/Hgi-j^Cd^rTe quantum 
wells (QWs) of a double quantum well (DQW) struc- 
ture. In the first section, we recall the k ■ p band struc- 
ture calculations for the confined states of an individual 
HgTe/Hgi_2:Cd2.Te QW. In the second section, we show 
how to calculate the tunneling matrix elements arising in 
a double quantum well configuration, starting from the 
knowledge of the envelope functions of a single QW and 
present the numerical values of the single-particle tun- 
neling amplitudes. 



1. Band structure model 

For the calculation of the band structure and wave 
functions of the single HgTe/Hgi-j^CdajTe QW, an 
envelope-function approximation^, based on an eight- 
band k ■ p Hamiltonian, is used. The total wave function 
is given as follows: 
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FIG. 4: (Color online) Band edge profile of a 

HgTe/Hgi_^Cd^Te QW of thickness d. The profiles of 
conduction and valence band are shown together with the 
spht-ofT band (with the split-off parameter A). 



layers: 



ui{r) = 


ire, +1/2) = 


= St 


U2{r) = 


ire, -1/2) = 


= Si 


U3{r) = 


|r8,+3/2> = 


= il/V2){X + tY)t 


Uiir) = 


|r8,+i/2> = 


= {l/V6)[{X + iY)i-2Zt] 


U5{r) = 


|r8,-i/2> = 


= -{l/V6)[{X-iY)t+2Zl 


ueir) = 


|r8,-3/2) = 


= -{l/V2){X-iY)i 


urir) = 


|r7,+i/2) = 


= {l/V3)[{X + iY)i+Zt] 
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|r7,-i/2> = 


= {l/V3)[iX~iY)t-Zi]. 
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The Hamiltonian for a QW grown along [001] direction 
is then given byJ^ 
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Here, the band structure parameters 71, 72, 73, /i = 
(73 — 72) /2, 7 = (73 + 72) /2, n and F describe remote 
band contributions; P is the Kane momentum matrix el- 
ement; Ec{z) and Ey{z) are the conduction and valence 
band edges, respectively; A is the spin-orbit splitting en- 
ergy. It should be noted, that the in-plane wave vector 
{kx,ky) is a good quantum number, but kz should be 
replaced by the operator k^ = —id/dz. 

The band structure parameters for HgTe and 
Hgi_a;Cdj;Te are considered as piecewise constant for 
each of the layers with an abrupt change at the interfaces, 
according to the band edge profile in Fig. U) Using the 
correct operator ordering in the Hamiltonian [Eas. (jA3[) 
and (jA4p ]. in accordance with the envelope- function ap- 
proach derived by Bur1j^, provides us with an unambigu- 
ous determination of the interface boundary conditions. 
A detailed description of the model, as well as the values 
of the band structure parameters for HgTe and CdTe, is 
given in Ref. Ilfl- 

Solving the eigenvalue problem Hi = Ei (with f an 
eight-component envelope function vector) we determine 
the envelope functions fn-M (z) and energy levels near 



icii = for a single QW. In particular, we are interested in 
the first conduction and valence subbands which for the 
QW width close to the critical value 6.3 nm are denoted 
by El and HI—. Other subbands are sufficiently apart, so 
that the system can be described, in a reasonable spectral 
ranged, by a 4-band effective Hamiltonian^, with the 
basis given by the El and HI subband eigenstates at 
k|| = i\El,+), \E1,-), |F1,+), \m,-)). Using the 
eight-band k-p model described above we calculate these 
states: 
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The summation index n usually runs over all eight ba- 
sis states (see Ea. (jA2[) ). But for k|| = only envelope 
functions with index n given in Eqs. (|A5| ) have nonzero 
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FIG. 5: (Color online) Envelope functions of the El and HI 
bands of a HgTe QW of thickness d = 5.7 nm, in panels (a) 
and (b), respectively, at ku = 0. 



values. In Fig. [5ja) and (b), we present the envelope 
functions for El and HI subbands at ky = 0. Note that 

/^k,ito(^) a'^d f8-k'n=oi^) giv'5 negligibly small contribu- 



El. 



cEl,+ 



tions to the total wave function (|(/7.ji-'^o('^)l/7-Jc'=o(^))l' 

\{fi',-oi^)\fil;=oim<oM). 



2. Tunneling matrix in a double quantum w^ell 

We have proposed an extension of the 4-band effective 
Hamiltonian^ for a HgTe/Hgi_^Cd:,Te QW to a DQW 
structure. Here, we present a realistic estimation of the 
parameters in the model equation Eq. [2l We consider a 
DQW geometry as described in Fig. El where we assume 
a symmetric structure for simplicity. Two HgTe QWs of 
individual thickness d are separated by a barrier of length 
t of Hgi_^Cd2:Te. When the barrier is sufficiently thin, 
the electronic states belonging to different layers share a 
finite overlap. The envelope function of the DQW can 
be described as a bonding and antibonding combination 
of the envelope functions of the individual f and b QWs. 
Each QW has two confined states in the spectral range 
of interest: the El and the HI subbands. Using the 
envelope functions for the El and the HI bands discussed 
in Section [XT] (see Egs. lASI and Fig.[S]), we can write the 
tunneling matrix between the states localized on f and b 
wells as 



hr = {X,+J\H\X,+,b) 
hr = {X,+J\n\X,-,b), 



(A6) 
(A7) 



with 



(r|X,±,6) = J2fn:k%iz-zb) ur^ir) (A8) 

n 

{r\X,±J) = Y.fn-M%i^-^f) Mr) (A9) 

n 

where zj — ~^^ and Zf, — i^ are the centers of the 
f and b QWs, and X refers to the El and the HI sub- 
bands. Ht is the tunneling matrix previously introduced 
in our model (in the form of Ht = hTVx), while hx is 
an additional Rashba-like term, mixing the two Kramer's 
block. H is the Hamiltonian of the DQW system, which, 
in the eight-band Kane model basis, has the same form 
as Eq. (|A3p and follows the band edge profiles of Fig. [51 
with its material-related parameters tuned to their HgTe 
values in the two intervals [—d — |, ~|] and [^,d + |], 
and to the Hgo.3Cdo.7Te values elsewhere. Note that we 
neglect the distortion of the band edges due to the elec- 
trostatic potential V^ which in our proposal contributes 
just a negligible correction with respect to the band off- 
sets between HgTe and Hgi_a;Cd2:Te . We stress that the 
results of any concrete calculation is obtained for a Cd 
content in the barriers of x = 0.7. The parameter x of- 
fers a way to tune the transparency of the barrier. In 
particular a smaller x leads to reduced band edge offsets 
between the QW and the barrier region (see Fig. ^. 
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FIG. 6: (Color online) Sketch of the band edge profiles of 
a HgTe/Hgi_2,Cda:Te DQW. The profiles of conduction and 
valence band are shown together with the split-off band (with 
the split-off parameter A), d is the thickness of a single HgTe 
QW and t is the barrier thickness of Hgi_a;Cd2,Te . 



Using the envelope functions of the El and HI sub- 
bands, calculated in Section [A 1[ we calculate the expec- 
tation values (integration over z) in Eq. (|A6|) in order to 
obtain hx- As expected from the band symmetry (see the 
related discussion in the note Il5l ) , hr has the following 
structure 



Abi 
a fc_ 



a k+ 
Affi 



(AlO) 



with 



Aei = 2 {El,+J\n\El,+,b) 

Ahi = 2(m,+,/|H|m,+,6) 

ak+ = 2 {El,+J\n\Hl,+,b) 



In Fig. [Tja), we show the dependence of Aei, Ahi 
and akx on the QW separation t for a fixed value of 
kx =0.1 nm^^ {ky =0). As expected, both A^i and 
a exponentially decay for increasing length of the tun- 
neling barrier. We note, also, that A^i is very small 
and for all practical purposes it can be taken as zero. 
The quantitative difference between Ahi and A^i is a 
direct consequence of the El band behaving as an interfa- 
cial confined state between HgTe and Hgi_j;Cda;Te , see 
Fig. [5l HI envelope function is, instead, mostly confined 
in the HgTe layer. 

In Fig. [7jb) and (c), we show akx and A^i as a func- 
tion of the in-plane wavevector kx for t = 7 nm. Param- 
eter Aei shows negligibly weak quadratic corrections, 
while a is strictly a constant. This justifies the treat- 
ment of the parameters a and A^i as constants in Eq.[2] 
Table H] provides the values of a and A^i for several val- 
ues of t, for a QW of width d — 5.7 nm. 

Here should be noted that BHZ model is applicable in 
a finite spectral range?-*^, this fact restricts the maximum 
value of the tunneling matrix elements for which the ap- 
proach we follow is valid. We performed calculations for 
d in a range between 5.7 nm and 6.6 nm, observing that 
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FIG. 7: (Color online) Tunneling terms of Eq. (|A10|) 
between the confined states of the b and f QWs of a 
HgTe/Hgi_a;Cd2,Te DQW structure. Tunneling terms are 
plotted as a function of t, for k^c = 0.1 nm^^ and ky — 0, 
in (a), while they are shown as a function of k^, for i = 7 nm, 
in (b) and (c). 



t (nm) 


2 


3 


4 


5 


6 


7 


Aisi (meV) 
a (meV nm) 


80 
66 


42 
32 


22 

18 


12 
10 


6 

5 


3.4 
3 



TABLE I: Tunneling parameters a and A_bi, calculated at 
k:c ~ ky = for several values of the barrier thickness t. d — 
5.7 nm, corresponding to a Dirac mass term M « 6.5 meV. 



the tunneling elements smoothly vary. In particular, for 
increasing d, Aei is slightly reduced (there is a varia- 
tion of about 10% between values corresponding to QW 
widths of 5.7 nm and 6.6 nm). Parameters a and Ahi 
are sensibly increasing with d, however in all the explored 
range A^i >> Ahi- 

The DQW geometry, breaking the single QW mirror 
symmetry, generates also Rashba-like tunneling terms de- 
scribed by Eq. (jA7p . which are not included in our model 
for clarity of the presentation. Our calculation shows 
that the only significant Rashba-like tunneling element is 
{El,+, flTilEl^—^b) = ffc-- However, in our proposal 
there is an energy potential shift V sa ±2Af between f and 
b QWs. This energy detuning between the El confined 
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FIG. 8: (Color online) Particular of the fist conduction and 
first valence bands of a DQW structure with t — 7 nm, d — 
5.7 nm, obtained from our model, and including the Rashba- 
like tunneling term in Eq. (|A7|) . Full lines and dashed lines 
refer to the two Kramer's blocks. 



t (nm) 



Ey (meV) 



6 



7 



-3.9 



-1.6 -0.85 



TABLE II: Bulk gap originated from the tunneling elements 
in the DQW structure with individual QW width d — 5.7 nm, 
corresponding to a Dirac mass term M ~ 6.5 meV. 



levels in the f and b QWs strongly hinders their cou- 
pling by the Rashba-like interaction. In practice, when 
we include the Rashba-like term in the model, using the 
realistic tunneling elements we have calculated, it results 
in a small (compared to the bandgap) tilting of the bulk 
dispersion curves (see Fig.|S|). Due to this, the dispersion 
curves are characterized by i?±(icjj) ^ i?±(— k||), with 
only the time reversal symmetry requirement E-^-(kn) = 
i?_(— k|j) still holding. We note that by turning off adi- 
abatically the Rashba-like term, we can connect the sys- 
tem with Rashba-like interactions to that previously pre- 
sented in the main body of this work without closing the 
bandgap. This proves that the topological properties of 
the system are not affected by the Rashba-like tunneling 
term in Eq. (ITf)) . 

In Fig. [SI we show the bulk dispersion curve of a 
DQW structure employing the tunneling parameters of 
t = 7 nm (see Table I), with M w 6.5 meV, corre- 
sponding to d = 5.7 nm, and with a potential shift 
T^ = — 15 meV (wc note that the electric field employed 
is about 1 mV/nm, which should be easily achievable in 
semiconductor heterostructures) . Full lines and dashed 
lines refer to the two Kramer's blocks. The Rashba-like 
term in Eq. (IA7|) is included, with a ~ 6 meV nm, but its 
effect, as previously explained, is negligible and consists 
in a slight tilting of the bands. 

A semiconducting gap of 0.85 meV arises, in agreement 
with the value obtained with the reduced model Eg [TUl 
using the formula 



E„ 



52 



{A^-AMB){D^-& 



(All) 



Table im shows the value of the bulk gap for several values 
oft. 



For i < 5 nm, the reduced model (and therefore 
Eq. lAllI) is no longer valid due to the large tunneling 
terms (Ht) compared to the individual QW mass term 



M, however the full model [Eq.(3)] still holds (as long 
as the BHZ model is a reasonable approximation of the 
band structure). 
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